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The random deposition model must be enriched to reflect the variety of surface rough- 
ness due to some material characteristics of the film growing by vacuum deposition or 
sputtering. The essence of the computer simulation in this case is to account for possible 
surface migration of atoms just after the deposition, in connection with binding energy 
between atoms (as the mechanism provoking the diffusion) and/or diffusion energy bar- 
rier. The interplay of these two factors leads to different morphologies of the growing 
surfaces from flat and smooth ones, to rough and spiky ones. In this paper we extended 
our earlier calculation by applying some extra diffusion barrier at the edges of terrace- 
like structures, known as Ehrlich-Schwoebel barrier. It is experimentally observed that 
atoms avoid descending when the terrace edge is approach and these barriers mimic this 
tendency. Results of our Monte Carlo computer simulations are discussed in terms of 
surface roughness, and compared with other model calculations and some experiments 
from literature. The power law of the surface roughness <r against film thickness t was 
confirmed. The nonzero minimum value of the growth exponent (3 near 0.2 was obtained 
which is due to the limited range of the surface diffusion and the Ehrlich-Schwoebel bar- 
rier. Observations for different diffusion range are also discussed. The results are also 
confronted with some deterministic growth models. 
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1. Introduction 

The surface growth phenomena are of interest for scientists because of thin films nu- 
merous industrial aspects as well as examples of that phenomena in crystal growth, 
biology etcJd'EI A great variety of models, approaches and approximations are de- 
voted to this problem, among them manv theories of solid film surfaces growth based 



on continuous as well as discrete models.ao The so called solid-on-solid (SOS) model 
is usually assumed. In this model no pores or voids or even overhangs are consid- 
ered and each particle sits on the top of other one. In SOS model the film surface 



1 



2 K. Malarz & A. Z. Maksymowicz 



morphology is described by single- valued function h(f, t) of surface height in planar 
coordinate r at time t. However, such information is usually too rich and not nec- 
essary. Instead, we characterize the surface in terms of some statistical parameters 
describing surface roughness, spatial correlations, anisotropy or growth dynamics. 
In Ref. a model for (l+l)-dimensional case is presented. We use this model for 
the (2+l)-dimensional case above the percolation limit when the substrate is cov- 
ered by at least one monolayer (ML) on average. We study surface morphology and 
time evolution of the surface width (spatial height-height auto-correlations). The 
results of simulations are compared with theoretical predictions as well as accessible 
experimental data. 

2. Model 

In the presented model particles simply fall down until they reach the top of a col- 
umn at randomly chosen position r = (x, y) on a two-dimensional lattice. Each par- 
ticle is represented by a unit- volume cube which can occupy only discrete position in 
the lattice. The simple cubic symmetry is assumed for simplicity of computations. 

To mimic molecular beam epitaxy (MBE) growth, one assumes that the dom- 
inant physical mechanism of the surface smoothing (or roughening) is surface dif- 
fusion. We assume that desorption from the surface can be neglected. We also 
simplify the situation to the SOS model where no overhangs or bulk vacancies are 
created. All these assumption are reasonably well satisfied in typical experimental 
situations. 

After initial contact to the surface, the particle is allowed to make several one 
step jumps to the nearest neighbors. We allow for Ldif jumps which defines the 
diffusion range. The probability of a particle to settle at a given site is given by the 
Boltzmann factor: p cx exp(—E/kT), where E is particle energy at the site, k is the 
Boltzmann constant and T denotes the temperature. There are three contributions 
to the energy E. The first one is the particle-particle binding energy J between the 
nearest neighbors. Then the diffusion barrier V enters if we claim that the particle 
jumps to a neighboring site. The third contribution S is an interaction energy 
between the top particle and underlying atom. The S'-term is an invention used in 
simulations to account for the edge effect. It was shownBCMa that atoms at the edge 
of a terrace are reluctant to step down. This situation was^imulated in theoretical 
study by an additional energy barrier for the edge atoms This barrier is known 
in the literature as the Ehrlich-Schwoebel (ES) barrier, although its origin is still 
not quite comprehended in terms of a physical mechanism directly responsible for 
such a barrier. Note, that in our model the ES-likc barrier does not make it difficult 
to hop down the step edge, but only to approach the step-edge site (similarly like 
in Ref. lid for instance) . 

Finally we get 

p(stay) oc exp — , (1) 
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Figure 1: A cross-section through the film growing on the substrate. Squares denote 
actually occupied sites in the lattice. The black squares are new arriving particles, 
(a) The gray squares are atoms which contribute to the ^-term. (b) The gray square 
is atom which contributes to J. (c) If the new (black) atom jumps to the nearest 
site, it must overcome the diffusion barrier V. The line of bricks indicates substrate 
and fat arrow shows direction of growth. 



and 

—(V + y J + y S) 

p(move left, right, forward, backward) oc exp — — — , (2) 

kT 

where the summation in y^ J is over nearest neighbors, and the sum yj S runs over 
nearest and next-nearest neighbors of the underlying layer, as shown in Fig. [t| 

So, the particle behavior after first contact to the surface depends on competition 
among values of J, S and V. (In calculations, we express them in kT units.) 
Negative values of J correspond to attractive forces between particles and, therefore, 
a tendency to make the surface smoother. Positive J describes repulsion. We assume 
negative values of S (S < 0). The value of the diffusion barrier is always positive 
(V > 0). With increasing V, the probability of diffusion to places marked by arrows 
in Fig. |l]c decreases. 

In Fig. H black squares are places of initial impact of particle to surface. The gray 
ones are possible places for particle incorporation after one step local relaxation. A 
height difference in positions (1), (2) and (3) in Fig. ||a (hi > h 2 > h 3 ) does not 
influence the preferences of the place of incorporation. However, for positive values 
of J, place (1) is more attractive than (2) and (3) for minimum energy principle is 
then equivalent to minimum contact to neighbors. And vice versa, for negative J 
sites (2) or (3) are preferred as offering the strongest bond. We must remember, 
however, that diffusion barrier V or the edge S-term may alter the preferences. 
In Fig. ||b both gray sites are equal energetically independently of the values of 
model control parameters. For non-zero value of S particles meet ES barrier when 
attempting position (3) in Fig. ||c. In all three cases in Fig. ^, the increase in 
diffusion barrier V promotes black positions (2). 



3. Results of Simulation 
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Figure 2: Squares denote actually occupied sites in the lattice. Gray sites denote 
possible places for incorporation to the surface, while the black square is the place 
of initial contact to the surface. The line of bricks indicates substrate and fat arrow 
shows direction of growth. 



Our results of simulations for the surface morphology may be compared with the 
scaling law. For a wide variety of discrete and continuum models of surface growth, 
the surface width a — wh ich d escribes surface roughness — is expected to follow the 
dynamical Family- VicsekE3£3 scaling law: 

<r(t,L) (xL a ■ f(t/L z ) (3a) 
... I x@ and z = a/(3 fori«l 

/(*) = <i , „ 1 ( 3b ) 

I 1 lor x >■ 1 

where a 2 = (h 2 ) — (h) 2 denotes the standard deviation of the surface height from 
the average surface height (h), L is the linear size of the system, t denotes time, 
the exponents a,0,z depend mainly on the system dimension, and (...) denotes 
configurational average over all sites r at time t. 

Scaling versus both time and space, with two characteristic scaling exponents: 
roughness exponent a and dynamic exponent z, results from the universal self-affine 
scaling. The surface profile and its properties are statistically invariant if the length 
in a direction parallel to the surface is scaled by a factor A and simultaneously 
the length in the perpendicular direction and the time t by factors A Q and A z , 
respectivelyEiB Such situation was also widely observed experimentally — see Ref. 
E-3 for review. 

3.1. Surface roughness 

The simulations were carried out mainly on 500 x 500 large square lattices with 
periodic boundary conditions. The total number of te n m onolayers were deposited. 

The scanning-tunneling microscopy (STM) studiesE£l show, that at higher tem- 
peratures particles can jump even over four lattice constants. However single or 
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double jumps are most frequent. Therefore, and for simplicity of computations we 
decided to assume single step relaxation, to the nearest neighbors only (£dif = 1). 
The results of simulations of surface growth for different control parameters J, S 
and V are presented in Tables 

3.1.1. Layer-by-layer growth 

For negative values of binding energy (J < 0), and for the diffusion barrier small 
enough to allow for the migration, a layer-by-layer (or Frank-van der Merve) growth 
takes place, when each layer is completed before aggregation of island on higher 
levels. 

In the limit J — > — oo a perfect wetting is expected for which a — 0, when all 
columns have the same heights (h). However, from the simulation we got that for 
smooth surface limit (presented in Fig. ^|) there are also the many sites with heights 
(h) ± 1 and some with (h) ± 2 (see histogram in Fig. |J). As result we get (3 w 0.22 
instead of j3 = 0. 

The deviation from ideally flat surface comes from limitation to one step only 
diffusion. The limit cr — ^ is for total energy minimum, while we seek local energy 
minimum only when employing local information on the site in question and its 
vicinity. Then increasing the diffusion range to Ldif is just Ldif times repetition of 
the same local search. This may either lead to the new final coordinates as far as 
Ldif steps away from the initially hit site, yet still before the energy minimum is 
reached, or the particle may oscillate to and fro around the hit site if it happens to 
be the local minimum. In each case, however, the final configuration is not the total 
energy minimum which would have been for the ideally flat surface. Thus for the 
finite ranges Ldif of diffusion, and for the applied algorithm of local search of energy 
minimum, we get in the J — > — oo limit er > 0. The dependence of cx(Ldif , J — * — oo) 
is discussed in chapter 3.1.5. 

Eq. (||) shows that surface properties depends strongly on system size L. Surface 
width a grows initially like t@ until a steady state and saturated roughness oc L a 
after some time r oc L z is reached. Determination of the roughness exponent a is 
much more difficult for large lattices because a huge number of layers must be 
deposited before the stationary state is reached, which requires enormous computer 
times. We were able to measure roughness exponent a only for not too large lattices 
(max. 64 x 64) and it seems that a « 0.78. 

3.1.2. Random deposition 

For J = S = (no driving force for particle migration) or for V — > +oo (when 
particles movements are frozen by the barrier), the surface growth corresponds to 
a simple random deposition model ( "sit and stay" ) . For J = S = we have 



p(stay) = p(left) = p(right) = p(forward) = p(backward) = 1/5, 
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Figure 3: Film profile for smooth surface (J — > — oo, (h)=10 ML). 
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Figure 4: Histogram of surface heights for smooth surface (J 
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Figure 5: Film profile for simple random deposition (J — 0, (h) — 10 ML). 
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Figure 6: Histogram of surface heights for random deposition (J = 0, (h)=10 ML) — 
Poisson distribution. 

while in the latter case V^oowe get 

p(stay) = 1 and p(move) = 0. 

The produced randomly surface profile is presented in Fig. 0. We recovered 
the well known result, of the random deposition model that the column heights 
follow the Poisson distribution (see histogram in Fig. That also implies a = \/t 
((3 = 1/2), the expected dependence of surface width a on time t (film thickness 
(h)). 

3.1.3. Growing of rough surfaces 

When jumps to the top of higher neighbor columnsare allowed, an unstable growth 
takes place and the surface width grows infinitely.Ej For J — > +oo the dependence 
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Figure 7: Film profile for rough, spiky surface (J — > +00, (ft)=10 ML). 
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Figure 8: Histogram of surface heights for rough surface (J — > +00, (h)=10 ML). 
N(h = 0) = 160000. 

of surface width against time becomes linear. Such rapid roughness growth (/3 w 1) 
was experimentally ohserved by atomic force microscopy (AFM) studies of plasma- 
etched silicon surfaces!!] and predicted by models of columnar growth.tS Note, that 
this time the "surface growth" means removing some atoms instead of depositing 
them. 

The surface profile for rough surface is presented in Fig. |?]. In such case there 
are mainly unoccupied sites (h = 0) while occupied sites have two or three times 
larger heights than the average surface height (h) (see Fig. ||). 

3.1.4. Surface profile dynamics 

The dynamical properties of film surface growth may be discussed in terms of growth 
exponent j3 (or dynamic exponent z) and roughness exponent a in Eq. (|). Nu- 
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Table 1: Surface width a [ML] for ten deposited layers, different values of diffusion 
barrier V and interaction energy J, with S = 0. 



V/J 


-5 


-1 


-0.5 





0.5 


1 


5 





0.87 


1.15 


1.52 


3.16 


7.34 


10.72 


14.17 


1 


0.88 


1.20 


1.58 


3.16 


7.59 


10.61 


13.90 


2 


0.88 


1.35 


1.80 


3.16 


6.94 


9.97 


13.70 


3 


0.89 


1.60 


2.17 


3.16 


5.45 


8.77 


13.44 


4 


0.92 


1.96 


2.58 


3.16 


4.16 


7.03 


13.13 


5 


0.96 


2.36 


2.88 


3.16 


3.52 


5.18 


12.71 


to 


1.39 


3.15 


3.16 


3.16 


3.17 


3.18 


11.29 


15 


1.98 


3.17 


3.16 


3.16 


3.16 


3.16 


9.29 


20 


2.61 


3.16 


3.17 


3.16 


3.16 


3.17 


5.57 


50 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.17 



Table 2: Surface width a [ML] for ten layers for different values of energies J and 
S, with no diffusion barrier V = 0. 



J/S 





-0.25 


-0.5 


-0.75 


-1 


-2 


-3 


-4 


-5 





3.16 


2.02 


1.55 


1.34 


1.23 


1.03 


0.97 


0.95 


0.95 


-0.25 


2.02 


1.49 


1.28 


1.16 


1.09 


0.97 


0.93 


0.92 


0.92 


-0.5 


1.52 


1.25 


1.13 


1.06 


1.02 


0.94 


0.92 


0.91 


0.90 


-0.75 


1.28 


1.13 


1.04 


0.99 


0.96 


0.91 


0.90 


0.90 


0.89 


-1 


1.15 


1.05 


0.99 


0.96 


0.94 


0.89 


0.89 


0.89 


0.89 


-2 


0.94 


0.92 


0.90 


0.89 


0.88 


0.87 


0.86 


0.87 


0.88 


-3 


0.90 


0.89 


0.88 


0.87 


0.86 


0.85 


0.86 


0.87 


0.87 


-4 


0.87 


0.87 


0.87 


0.86 


0.86 


0.85 


0.86 


0.86 


0.86 


-5 


0.88 


0.87 


0.87 


0.86 


0.87 


0.85 


0.85 


0.86 


0.86 



Table 3: Surface width a [ML] for ten layers for different values of diffusion barrier 
V and energy S, with J = 0. 



v/s 





-0.25 


-0.5 


-0.75 


-1 


-2 


-3 


-4 


-5 





3.16 


2.02 


1.56 


1.34 


1.23 


1.03 


0.97 


0.95 


0.95 


1 


3.15 


2.07 


1.63 


1.40 


1.28 


1.06 


0.98 


0.96 


0.95 


5 


3.16 


3.04 


2.86 


2.62 


2.34 


1.62 


1.36 


1.22 


1.11 


10 


3.16 


3.16 


3.17 


3.15 


3.15 


2.82 


2.10 


1.76 


1.52 


50 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 
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Table 4: Surface width a [ML] for ten layers for different one-step diffusion ranges 
L dif and binding energy J (V = and 5 = 0). In all cases <r(L dl{ ) decreases with 
increasing Ldif until some saturated roughness <T sat is reached. 



Ldif/ J 


-20 


-5 


-1 


-0.75 


-0.5 


-0.25 


-0.1 


-0.05 








3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


3.16 


1 


0.87 


0.87 


1.15 


1.28 


1.52 


2.02 


2.61 


2.86 


3.16 


2 


0.82 


0.81 


1.03 


1.17 


1.42 


1.98 


2.61 


2.87 


3.16 


5 


0.80 


0.80 


0.94 


1.06 


1.29 


1.85 


2.54 


2.83 


3.16 


50 


0.68 


0.69 


0.93 


1.04 


1.26 


1.83 


2.52 


2.82 


3.16 


500 


0.47 


0.50 


1.00 


1.10 


1.30 


1.83 


2.53 


2.82 


3.16 


5000 


0.36 


0.49 


1.04 


1.11 


1.30 


1.83 


2.52 


2.84 


3.16 



merically we found that growth exponent (3 depends roughly only on the sign of 
interaction energy J (when S = V = 0): 

{0.22 for J — > -oo 
1/2 forJ->0 
1 for J — * +oo 

independently on lattice size L. 

Before completing the first monolayer a deviation from the true value of f3 to- 
wards 1/2 could be observed as an artifact of the initially perfectly flat substrate. 
Our results axe consistent with some deterministic SOS models studied by Kotrla 
and Smilauenlj (where 0.25 < /? < 1), Amar and FamjlyEJ SOS models for bcc 
structure (/3 ~ 0.24), and He et al experimental estimateEil (/3 ss 0.22). 

3.1.5. Diffusion range 

So far presented results were produced for diffusion limited only to one inter-atomic 
separation (Ldif = 1)- It was earlier suggested by FamilyEl that growth mode is 
unchanged independently of diffusion range yet it is finite. Here we check how diffu- 
sion range Ldif influence on surface morphology for the case adequate to the Family 
model, it means J — > — oo. We found that an increasing diffusion range decreases 
surface roughness (see Table |]) until some saturation level <r S at is reached. The 
saturation roughness level depends on binding energy J. The random deposition 
case is recovered for Ldif = and/or J = 0. 

3.1.6. Substrate temperature 

Table || shows how the substrate temperature alters the surface roughness. The 
values of model control parameters were taken after Ref. 11 where simulation based 
on Arrhenius dynamics for homoepitaxy on Pt(lll) substrate were carried out. Here 
we identify our energetic diffusion barrier V as the activation for diffusion energy 
(V = 0.75 eV), the interaction energy J as activation for breaking in-plane bonds 
between Pt atoms energy (J = —0.15 eV), and activation for breaking bonds with 
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Table 5: Surface width a after deposition of 20 layers for different substrate tempera- 
tures T for homo-epitaxy growth on Pt(lll) substrate: V = 0.75 eV, J = —0.18 eV, 

S = -0.15 eV. 

T [K] 300 600 900 1200 1500 
a [ML] 2.04 1.87 1.74 1.60 1.49 



substrate energy [S = —0.18 eV) as our S. The increase in temperature T makes 
particles more mobile and then they may escape bonds easier, and across energetic 
barriers, which results in a decrease of surface roughness a. Such a situation is 



typical and well verified both experimental 
relation was reported for instance in Ref. 

4. Conclusions 



and theoretically^ However, no simple 



We obtain from the simulation 0.22 < (3 < 1 and that lower limit (3 w 0.22 is caused 
by the limited diffusions which pushes (3 from (3 = (perfect wetting) towards 
(3=1/2 (random deposition). 

We confirmed results from simulations based on deterministic rules of particles 
relaxationllj that allowing or preventing particles to climb to higher levels is essential 
for extracting the growth and roughness properties. When such jumps are forbidden, 
then a grows accordingly to Family- Vicsek law (||). In contrast, when jumps to the 
top of higher neighbor columns are allowed, unstable growth takes place and surface 
width grows infinitely. 

Indeed, from the time evolution of surface profiles we confirm that for negative 
values of J (when particles attract each other and as a result of that maximize the 
number of bonds like in the Wolf- Villain modelE3 — or minimize their heights like in 
the Family modeltj) a stationary state is observed (J — * — oo,a » 0.78, (3 w 0.22). 
The saturation of roughness (lim^oo a = <7oo) as a function of time (total coverage) 
was reported also experimentally for instance in Liu et al worked for gold films 
growth on mica substrate under different experimental conditions. 

In the other limit, when particles repel each other and jump to higher levels 
(against gravitational force) the surface width growths infinitely and no station- 
ary state is observed; this could be considered as infinite values of the roughness 
exponent a ( J — * +oo,a — > +00,(3 = 1). 

For finite ranges of diffusion, there exists some saturated roughness which does 
not allow for perfect wetting observed only in limiting case: 

lim lim cr sat = 0. 

J — > — OO -Ldif — >oo 
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